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We theoretically demonstrate the possibility to observe the macroscopic Zeno effect for nonlinear 
waveguides with a localized dissipation. We show the existence of stable stationary flows, which are 
balanced by the losses in the dissipative domain. The macroscopic Zeno effect manifests itself in the 
non-monotonic dependence of the stationary flow on the strength of the dissipation. In particular, 
we highlight the importance of the parameters of the dissipation to observe the phenomenon. Our 
results are applicable to a large variety of systems, including condensates of atoms or quasi-particles 
and optical waveguides. 
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Since the pioneering work of Khalfin concerning the 
non-exponential decay of unstable atoms [l| the relation 
between the decay rate and the measurement process was 
in the focus of many studies. One of the fundamental re- 
sults of the theory, termed after the seminal paper Q, 
as quantum Zeno effect, consists in slowing down the dy- 
namics of a quantum system subjected to frequent mea- 
surements or to a strong coupling to another quantum 
system. This phenomenon was demonstrated in a rigor- 
ous mathematical framework in Q and received its fur- 
ther refinements and extension in subsequent studies Q . 
Experimentally, the quantum Zeno effect has been con- 
firmed for single ions ultracold atoms in accelerated 
optical lattices [5] , atomic spin motion controlled by cir- 
cularly polarized light [6j], an externally driven mixture 
of two hyperfine states of neutral atoms [7] , photons in a 
cavity [8] and the production of cold molecular gase s in 
an optical lattice [9] . It has also been predicted jiol [Til ] 
that the tunneling dynamics of particles in a double- well 
potential can be slowed down if the particles are removed 
from one of the wells. Qualitatively similar results for 
the suppression of atom losses in an open Bose-Hubbard 
chain were reported in [l2[. In the limit of an infinitely 
strong measurement of particles in a given spatial domain 
it has been shown that the system is projected onto a 
unitary dynamics in the loss- free domain [13j . 

The Zeno effect is sometimes also understood in more 
general terms as the effects of changing a decay law de- 
pending on the frequency of measurements [14[. Apply- 
ing this definition to a macroscopic quantum system, like 
a gas of condensed bosonic atoms and, taking into ac- 
count that in the macroscopic dynamics the frequency of 
measurement can be interpreted as the strength of the 
induced dissipation the effect of the measurement 
on the decay of the quantum system can be viewed as 
the effect of dissipation on the macroscopic characteris- 



tics of the system. Here we assume this interpretation of 
the phenomenon and address the questions how the ap- 
pearance of localized losses in a waveguide is connected 
to the appearance of Zeno-like dynamics. In order to 
emphasize the distinction of the latter statement of the 
problem with respect to already standard and widely ac- 
cepted notion of quantum Zeno effect, below we refer 
to macroscopic Zeno effect (MZE) bearing in mind its 
meanfield manifestation. 

Losses are ubiquitous for real quantum systems due to 
the coupling to an environment. Very often, the loss pro- 
cesses are also spatially localized. They can be either ex- 
ternally engineered, e.g. with the tip of a scanning probe 
microscope, a local probe in a quantum gas, an absorbing 
spatial domain, or they can be intrinsically present in the 
form of defects and impurities. One can therefore expect, 
that the MZE can manifest itself in a wide class of phys- 
ical systems, including exciton-polaritons fl5j , magnon 
gases [l6j| 5 surface plasmons [ItJ, and optics of nonlinear 
Kerr media [l8| . 

We study one-dimensional nonlinear waveguides gov- 
erned by the nonlinear Schrodinger equation 



i^t = -^ xx + g\^\ 2 ^ - i-y(x)V, 



(1) 



where g the nonlinearity parameter and the local loss 
processes are modelled by ij(x) (for a review on the ap- 
plication of complex potentials see e.g. [19|). Since the 
localized dissipation is applied to a homogeneous con- 
densate, i.e. it breaks the translational invariance of the 
system, it can be referred to as a dissipative defect [ioj . 
We are interested in stationary flows, which correspond 
to a situation when an incoming flux of a particles from 
both ends of the waveguide is exactly balanced by the 
losses in the dissipative domain. Notice that in such a 
statement counting of lost particles is replaced by com- 
puting the number of particles that must be loaded into 
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FIG. 1: Possible experimental scenarios to observe the MZE: 
nonlinear optical waveguide (upper left), a magnon waveg- 
uide (upper right), a plasmonic nanostructure (lower left), an 
atomic BEC in a waveguide and two reservoirs (lower right). 



the system in order to compensate the losses. 

In Fig. [I] we show several examples of how such a sce- 
nario can be realized experimentally in different physical 
systems. While our approach is applicable to a large va- 
riety of physical situations, the systems we have in mind 
are those described in [2CH23j . i.e. an atomic BEC sub- 
jected to removal of atoms. In this last case, the time 
and the coordinate are respectively measured in units of 
(2uj±)~ 1 and a±/2, where a± and uj± are the transverse 
linear oscillator length and the frequency of the trans- 
verse trap, while g = a s no, being a s the scattering length 
and no the unperturbed linear density of the condensate. 
In the following we consider only the case g > which de- 
scribes repulsive inter-atomic interactions (or defocusing 
Kerr media in optical applications). 

The dissipation is described by the nonnegative local- 
ized function 7(2;), which is characterized by two control 
parameters: its amplitude To and characteristic width t. 
It is convenient to set j(x) = Tof(x/£) where f(x) is a 
known smooth function such that max^ \f(x)\=f(0)~l 
and maxj \f x (x)\ ~ 1. Then To and £ are proportional 
to the intensity of the defect: j(x)dx oc Tq£. We 
will also assume the most typical experimental situation 
where j(x) is an even function: j(x) = j(—x) with only 
one maximum at x = 0. Having in mind the experiments 
of Ref. [2ll-[23| one can estimate that To ~ Icrion/(eou±) 
where / is the current of the electron beam, eo is the elec- 
tric charge of the electron, and <7^ on is the total ionization 
cross section. 

The stationary flows are sought in the form ^ s t(t, x) = 
p(x) exp [i v(s)ds — ipt\ , where v(x) is the superfluid 
velocity, p is the chemical potential, and p 2 (x) = n(x) is 
the density. Substituting ^ st (t,x) into Eq. (pp) we obtain 
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0, j x + ~f(x)p 2 = (2) 



with j(x) = v(x)n(x) being the superfluid current. We 



are interested in solutions of Eqs. ((2]) with a constant 
density at infinity: lim^^^ \p(x)\ = p^. Then p = 
jloP^ + 9ploi where joo = =p lim^ioo j(x) is a positive 
constant. For any stationary flow, the loss of particles 
in the defect has to be balanced by the incoming current 
joo. The main objective of the present study is to show 
the existence of such stationary flows and to explore the 
dependence of the current joo on the parameters of the 
defect. 

First, we consider an example that allows for an 
exact solution, extending the result of [20|. We 
assume a dissipative defect of the particular form 
j(x) = 3r sech 2 (x/£). Then it is a straightforward to 
show that p{x) = tanh(x/£) and j(x) = —To£ta > nh 3 (x/£) 
are solutions of the system (j2j) provided that £ 2 (g + 
Tq£ 2 ) — 2. Thus the incoming flux is linearly propor- 
tional to the intensity of the dissipative defect: joo = To^, 
and in order to obtain a stationary solution, increase of 
the strength of the dissipation must be compensated with 
an increase of the incoming flux. In other words, if the 
incoming flux of particles is increased, the excess parti- 
cles can only be removed by a stronger defect. While 
this result is quite intuitive, we show below that it does 
not hold in general. In particular, we will show that for 
appropriate parameters, an increasing flux can be com- 
pensated by a weaker defect. 

We now focus on a dissipation with finite support: 
j(x) = if \x\ > £. This form of the dissipative term 
models, in particular, the electronic beam used in [2lll22|. 
In order to decrease numerical errors we choose j(x) 
to be smooth at the edges of the dissipative domain: 
7(0?) = r (1 -x 2 /£ 2 ) 2 if \x\ <£. 

Since Eq. ((2]) is not integrable - unlike its conservative 
counterpart where j(x) = - it is convenient to treat 
To as a parameter that increases departing from zero. 
Experimentally, this would correspond to an adiabatic 
increase of the defect intensity. For To = one recovers 
two well-known solutions: a constant density p{x) = p^ 
and a dark soliton p{x) = Poo tanh(^/ g/2 poox) [j(x) = 
for both solutions]. When the defect is adiabatically 
switched on, the constant density and the dark soliton 
give origin to two branches of solutions. The branch bi- 
furcating from the constant density consists of symmet- 
ric flows, for which the relation p{x) = p(—x) holds. The 
flows that branch off from the dark soliton are antisym- 
metric, p(x) = —p{—x). From the second of Eqs. (j2j) it 
follows that both the symmetric and antisymmetric flows 
possess odd currents: j(x) — —j(—x). 

Considering behavior of the solutions in the vicinity 
of x = 0, for the symmetric flows we obtain p xx (0) = 
P(°)[#0 2 (°) - Plo) ~ fooP^Y Thus > employing a physi- 
cally obvious condition p^ > p(0) > 0, we find p xx (0) < 
0. We therefore arrive at the counterintuitive conclusion 
that for the symmetric flows the atomic density n(x) has 
a local maximum in the point of maximal dissipation. 

On the other hand, for x — >• 00 both symmetric and 
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antisymmetric flows behave as 
Poo - p{x) oc e"^, where A = 2gp 2 00 - 4j 2 p- 4 . (3) 

Thus, for a given density p^ there exists an upper bound 
for the maximal current j™ ax = ^fgp^/V^ above which 
no stationary flow can exist. 

In Fig. [2] (a) [Fig. [2] (b)] we show the density profiles 
n{x) of symmetric [antisymmetric] flows for different val- 
ues of the dissipation strength r . We observe that for 
the symmetric flows, the density possesses two deep lo- 
cal minima. For a weak dissipation (e.g. Tq = 0.01) 
the minima are situated outside the dissipative domain, 
which allows to compute the exact value of the density in 
these minima: n m i n = 2j^ p^ 4 /g. As the strength of the 
dissipation grows, the minima move from ±oo towards 
the center and eventually enter the dissipative domain. 
For antisymmetric flows the dependence of the density on 
r is much weaker pronounced [the curves for different r 
are hardly distinguishable on the scale of the Fig. [2] (b)]. 

Now our goal is to study the dependence joo vs r . The 
typical results for the symmetric [antisymmetric] flows 
are illustrated in Fig. [2] (c) [Fig. [2] (d)] . When the non- 
linearity is strong enough (g = 1), for both types of the 
flows one can clearly see a global maximum of joo(ro)- 
When To exceeds the value corresponding to this max- 
imum, the required current joo decreases. This is the 
manifestation of the MZE. 

In order to observe the MZE in an experiment, it is 
important that the solutions are stable. To explore the 
stability of the flows, we substitute x) = ty st (t,x) + 
e" i/i£ [a+(x)e xt + a^_(x)e x * t ] into Eq. (pQ), linearize it 
with respect to a±(x), and solve the obtained linear 
eigenvalue problem. Instability occurs if there exists an 
eigenvalue A with positive real part A r . The results of the 
analysis are shown in Fig. [2] (c)-(e). For the symmetric 
flows presented in Fig.[2](c), small values of To do not al- 
low for stable stationary flows (see the dotted fragments 
of the lines). For instance, the symmetric flow with two 
local minima situated outside of the dissipative domain 
[shown in Fig. [2] (a) and corresponding to To = 0.01] is 
unstable. For larger values of To the symmetric flows be- 
come stable. We have observed that the two minima of a 
stable symmetric flow are always located in the dissipa- 
tive domain. All antisymmetric flows shown in Fig. [2] (d) 
are stable. Most importantly, the parameter range, in 
which the MZE is observed, has only stable flows. From 
experimental point of view, the existence of stable sym- 
metric flows is very appealing, as symmetric flows arise 
from the overall (symmetric) ground state of the system. 

We now elaborate on the role of the size of the defect. 
In Fig. [2] (e), we show the currents vs i for a fixed To- 
As the defect becomes wider, one might expect a mono- 
tonic increase of joo- As can be seen from the graph, 
this behavior is indeed encountered, however only in av- 
erage. Locally it is superimposed by a resonance-like 
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FIG. 2: (a): Density distributions n(x) for symmetric flows 
for g — 1 and £ — 4 and different values of IV Solid line: To = 

0. 01, dashed line: r = 1, dotted line: T = 10. (b): Density 
distributions n(x) for antisymmetric flows for g = 1 and £ = 

1. Solid line: To = 0.1, dashed line: ITo = 1, dotted line: 
To = 10. (c) and (d): Current vs strength of the dissipation 
for symmetric flows (with £ = 4) and for antisymmetric flows 
(with £ = 1) obtained for g = 0.1 and g = 1; stable (unstable) 
flows correspond to the solid (dotted) fragments of the curves; 
(e) : Currents and instability increments vs width of the defect 
for symmetric [(s)] and antisymmetric [(a)] flows for g = 1 and 
To = 1. In all panels poo — 1. 



structure which tends to enhance the current for certain 
defect sizes. We also find that stable solutions only ap- 
pear when the joo(£) is a growing function and that (up 
to a certain degree of accuracy) the domains of stabil- 
ity of the symmetric flows coincide with the domains of 
instability of the antisymmetric flows and vice versa. 

So far, we have encountered two different situations: 
in the case of the sech 2 -shaped dissipation, the branch 
of stationary flows does not show the MZE. In the case 
of the dissipation with finite support, the MZE has var- 
ious manifestations. The difference can be explained by 
looking at the asymptotic behavior of the corresponding 
flows. In the case of the sech 2 -shaped dissipation, the 
asymptotical behavior at x — >> oo of the density is com- 
pletely determined by the characteristic width t of the 
defect: p^ — p(x) ~ 2e~ 2x ^ , while the flows supported 
by the dissipation with finite support behave according 
to Eqs. (j3]). Moreover, Eqs. (j3]) imply that there exists 
the maximal possible current j™ ax . It appears that the 
presence of such a threshold is a signature of the MZE. 
No such threshold exists for the stationary flows with the 



sech 2 -shaped dissipation: the current joo can be arbitrar- 
ily large and the MZE is not found. 

However, the sech 2 -shaped dissipation still does not 
forbid the MZE in principle, since flows obeying Eqs. (|3j) 
can also be found in this case. Let us revisit the dissi- 
pation of the form j(x) = 3r sech 2 (x/^). Substituting 
p(x) = poo — pi(x), where pi(x) = o(l) as x — >> oo, into 
Eqs. (|2|) and neglecting the terms of smaller order, one 
observes that for x ^> 1 the function p\{x) is described 
by the equation p ljXX - Ap 1 = 12j oo T £p^e~ 2x/£ . For 
pi(x) to obey the asymptotics (j3]) the two conditions 
must be fulfilled: (i) A > - gives the maximal cur- 
rent: joo < j™ ax ; (h) a/A < 2/£ - yields the minimal 
possible current: > j™ in = poo^gp^/2 - £~ 2 (if the 
expression under the radical is negative, then j™ in = 0). 
As £ grows, j™ in approaches j™ ax . Hence, the range of 
currents allowing for the solutions which obey Eqs. (j3j) 
decreases. This leads us to the conjecture that rapidly de- 
caying dissipation is favorable for the observation of the 
macroscopic Zeno effect. In particular, the defects decay- 
ing faster than exponentially are more likely to display 
the MZE than ones obeying exponential decaying. 

We close this Letter with a discussion of possible ex- 
perimental observation of the found MZE. The incoming 
flux of particles has to be generated at both ends of the 
waveguide. This can be achieved by controlled pump- 
ing terms in the case of quasiparticles or by reservoirs 
in the case of real particles. For light propagating in a 
non-linear waveguide, such boundary conditions appear 
rather naturally. But even for a finite system with no 
reservoir, one can speculate that a quasi stationary state 
is established on intermediate time scales in a transient 
regime: if the defect is switched on in a finite system that 
is initially in its ground state, the stationary flow will de- 
velop out of symmetric initial conditions and will retain 
its symmetry with increasing dissipation. With time, a 
flow of particles towards the defect is created which can 
mimic the boundary conditions, applied in the preceding 
discussion. The condition of having a defect that drops 
faster than exponentially can be easily realized in most 
of the experimental implementations suggested above. 

We now support this reasoning by illustrating the gen- 
eration of stationary flows through direct integration of 
Eq. (pQ) on a finite domain subject to the boundary condi- 
tions ±L) = poo e ~ %llt (here L is the half- width of the 
computational domain). In Fig. [3] we show the temporal 
evolution of the atomic density for three different widths 
of the defect. For all the shown evolutions, the initial 
density is taken to be constant and the chosen bound- 
ary conditions fix the density and the chemical potential 
at the edges of the computational domain. Figure [3] (a) 
shows the evolution for a set of parameters, where a sta- 
ble symmetric stationary flow exists. After initial de- 
crease of the density at the location of the defect, the 
system achieves the stationary flow. In the vicinity of the 
origin one clearly observes the two local minima residing 




FIG. 3: Evolution of the density |^(£,x)| 2 starting from the 
initial data \£(0,x) = 1. For all the shown panels poo = 
g = To = 1. (a): The width of the defect is £ = 4. The 
generation of the symmetric stationary flow occurs, (b) and 
(c): The width of the defect is £ — 2 and £ = 6 respectively. 
The symmetric flows are unstable, and therefore no stationary 
flow is established. 



inside the defect. Indeed, the density approaches eventu- 
ally the stationary flow profile not only in the vicinity of 
the origin but in the entire computational domain. How- 
ever, if the symmetric flow is not expected to be stable 
[Fig. [3](b) and Fig. E9(c)], then the density profile shows 
ongoing distortions in the vicinity of the origin and even- 
tually loses the symmetry. For the cases shown in panels 
(b) and (c), the density tends to approach a profile corre- 
sponding to a stable antisymmetric flow which exists for 
the chosen values of £ and To [recall that the domains of 
the stability of the symmetric and antisymmetric flows 
alternate as shown in Fig. [2] (e)]. However, in the cases 
(b) and (c) truly stationary antisymmetric flows are not 
established because the chosen boundary conditions can 
support only symmetric flows, which are unstable for the 
chosen parameters. 

In summary, we have analyzed a non-linear waveguide 
with a localized dissipative defect. We found evidence for 
the appearance of the macroscopic Zeno effect for spe- 
cific boundary conditions and have analyzed to role of 
the interaction and the defect size. The observed macro- 
scopic Zeno effect is intrinsically related to the existence 
of the stable stationary flows and is expressed through 
nonmonotonic dependence of the intensity of the dissi- 
pation on the density of incoming currents. The proofed 
existence of stable solutions for symmetric flows is very 
important, as these solutions correspond to many natural 
experimental situations. 
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